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Electric double layer supercapacitors are promising devices for high-power energy storage based 
on the reversible absorption of ions into porous, conducting electrodes. Graphene is a particularly 
good candidate for the electrode material in supercapacitors due to its high conductivity and large 
surface area. In this paper we consider supercapacitor electrodes made from a stack of graphene 
sheets with randomly-inserted "spacer" molecules. We show that the large volumetric capacitances 
C > 100 F/cm^ observed experimentally can be understood as a result of collective intercalation of 
ions into the graphene stack and the accompanying nonlinear screening by graphene electrons that 
renormalizes the charge of the ion clusters. 



I. INTRODUCTION 

Electric double layer supercapacitors (SCs) are a 
promising class of efficient, long-lasting, high-power elec- 
trical energy storage devices based on reversible adsorp- 
tion of ions onto the surface of a porous electrode. The 
most common SC devices use activated carbons as elec- 
trode materials because of their high specific surface area 
and moderate cost. Recent improvements in electrode 
materials have greatly improved the energy storage ca- 
pacity of SCs, so that SCs are fast becoming viable can- 
didates to replace conventional batteries in a number of 
energy storage applications [H |2j . 

In this paper we study SC devices made from graphene 
electrodes, which have attracted considerable attention in 
recent years [SHS]. In such devices individual graphene 
layers are stacked to form an electrode and placed in con- 
tact with a reservoir of ionic liquid (or some other con- 
centrated electrolyte solution), as shown schematically in 
Fig. [ij A functional SC device contains two such elec- 
trodes with a voltage applied between them, usually with 
a porous separator between them that isolates the elec- 
trodes from each other electronically while allowing ions 
to flow from the reservoir to either electrode. 

Experiments on graphene-based SCs have demon- 
strated [21 El [5] volumetric capacitances in excess of 
lOOF/cm'^. These large capacitances are usually ex- 
plained within the simple picture in which the capacitor 
energy is stored in the electric double layer (EDL) that 
forms when ions are adsorbed onto the electrode's active 
surface. The volumetric capacitance of the electrode is 
then written as 

C = A^^Yii, , (1) 

where "^edl is the EDL capacitance per unit area and 
A is the surface area per unit volume of the electrode. 
In order to make a conservative estimate of the value 
of "^^EDL necessary to explain the large experimental val- 
ues of capacitance, one can imagine that an EDL forms 
on both sides of every graphene sheet and that adjacent 
graphene sheets have just enough separation to allow ions 




FIG. 1. (Color online) A schematic depiction of a graphene 
supercapacitor electrode. Graphene sheets are arranged in a 
stack that is placed in contact with a reservoir of ionic liquid. 
As a voltage is applied between this electrode (which we take 
to be the negative electrode) and the opposite electrode (not 
shown), cations (red circles with -|-'s) are driven to intercalate 
between graphene layers to neutralize the electrode's negative 
electronic charge. 

to fit sterically between them. In this case for an ion di- 
ameter a = 1 nm one finds that the value 100 F/cm^ im- 
plies an EDL capacitance of at least "^^edl ~ 7 //F/cm^. 
This value of "^edl is similar to the reported EDL capac- 
itances of a metal/ionic liquid interface [TUIIII]. Thus, to 
explain the large experimental value of C one should ap- 
parently imagine that each graphene sheet provides two 
independent EDLs whose capacitance is as large as that 
of a free metal/ionic liquid interface. 

Such large capacitance of graphene SCs is difficult to 
understand if one recalls that graphene is not a metal but 
a semi-metal. Graphene has a relatively small thermo- 
dynamic density of states (TDOS) v = drie/d^e, where 
TT-e is the two-dimensional (2D) concentration of electrons 
and /le is their chemical potential, which is tuned by the 
capacitor voltage. It is well known P^HTS] that within 
the mean-field theory this finite TDOS modifies the EDL 
differential capacitance as follows: 

"^EDL = '^EDL,oo + "^(T^ ' % = ■ (2) 

That is, one effectively has a "quantum capacitance" 
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"ioq that adds in series with the ideal EDL capacitance 
"^EDL.oo which is achievable only at the surface of a hy- 
pothetical perfect metal with infinite TDOS. Hence, the 
quantum capacitance imposes an upper bound on the 
volumetric capacitance of the system. 

This statement has profound implications for graphene 
where the TDOS has the form 
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where ksT is the thermal energy and v — 1.0 x 10^ m/s is 
the Fermi velocity. The corresponding quantum capaci- 
tance has a deep minimum at the neutrality point /ig = 
equal to 
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ksT = 0.83 ^F/cm^ 
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at T = 295 K. Therefore, the mean-field theory of Eq. Q 
predicts i) a minimum value of C that is a factor of eight 
smaller than what is observed and ii) a strongly varying 
/7-shaped voltage dependence of the capacitance, which 
is at odds with the roughly constant measured C{V). 

These same inconsistencies have been discussed in 
Ref.fTSl where an EDL capacitance "^edl = 7-10 ^F/cm^ 
was measured at the interface of a single graphene sheet 
and ionic liquid. It was suggested that graphene in this 
experiment was subject to randomly positioned charged 
impurities, causing smearing of the ideal TDOS and an 
increased '^q,min- By analogy, one can argue that multi- 
layered graphene SC devices [3HS]j which are doubtless 
heavily disordered, also have higher v than Eq. ([3| pre- 
dicts [T7] . 

In this paper we explore a more intriguing possibil- 
ity, namely that large enhancement of capacitance can 
result from a breakdown of the mean-field theory. We 
propose a theoretical model in which ions enter the elec- 
trode cooperatively as dense clusters, in analogy to stag- 
ing in intercalated graphite compounds [HI HH] • The ef- 
fective attraction between the like-charged ions is medi- 
ated by elastic stresses induced in the graphene stack fSU!- 
I24j . The high charge concentration inside such ion clus- 
ters activates a strong nonlinear screening in layered 
graphene p5H?7] that renormalizes the ion charge. We 
show that as a result the volumetric capacitance is greatly 
enhanced above the mean-field value. 

Our general approach to calculating the capacitance of 
a given electrode (say, the negative electrode) is as fol- 
lows. We first compute the total free energy F{Nj^.) per 
unit volume associated with the lowest energy configu- 
ration of cations (and the neutralizing concentration 
of electrons) per unit volume in the electrode. The 
value of the electronic charge Q = eN^ per unit volume 
of the electrode is that which minimizes the system's to- 
tal free energy density F — QV, where the term —QV 
represents the work done by the voltage source. Using 
the equilibrium condition d{F — QV)/dQ = gives 
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The differential capacitance per unit volume of the elec- 
trode, C — {dV/dQ)~^ , can therefore be written 
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In this way, the capacitance is closely related to the 
charge compressibility {N^d'^ F / dN^)~^ of the system: 
large capacitance implies high compressibility. The ca- 
pacitance can be expressed as a function of voltage, 
C{V), by combining Eqs. ([5| and These relations 
give the capacitance of a single electrode (here, the neg- 
ative electrode); the total capacitance of the SC device 
is the series sum of the anode and cathode capacitances. 
In this paper we concern ourselves with the capacitance 
of a single electrode only. For the majority of this paper 
we neglect entropic effects, so that in Eqs. (5]) and ^ 
the free energy F can be replaced by the total energy U. 

The remainder of this paper is organized as follows. 
Secti on [TT] defines the model system to be studied. In 
Sec. |III| we discuss the capacitance of the graphene 
electrode in the mean-field approximation. Section |IV| 
overviews the physics of staging and nonlinear screening 
in graphite intercalation compounds and calculates their 
capacitance. In Sec.|V]the mean-field and staging results 
for capacitance are united and the main results of this 
paper are derived. We show that the capacitance of the 
SC everywhere exceeds the mean-field quantum capaci- 
tance. Finally, in Sec. |VI| we provide additional discus- 
sion of the large-voltage regime, at which the mean-field 
quantum capacitance manifests itself in a somewhat en- 
hanced form, and we provide some concluding remarks. 



II. GRAPHENE-BASED SUPERCAPACITOR 
ELECTRODES 

Over the past few years a number of studies have 
confirmed the large capacitance per unit volume of SC 
devices in which the electrode is made from a stack 
of graphene sheets [SHSl |28j, as depicted in Fig. [T] 
Most commonly, in such devices the graphene sheets are 
ordered in the transverse direction, but have random 
translations or rotations between adjacent sheets ("tur- 
bostratic disorder"), so that on average the electronic 
dispersion relation of each graphene sheet is unaltered 
by inter-layer coupling. Below we refer to this arrange- 
ment as a "graphene stack" (GS). In practice, the layered 
ordering of graphene sheets within the GS electrode is 
preserved only over some mesoscopic length scale, usu- 
ally on the order of tens of nanometers [SI [55] . In the 
remainder of this paper we make the assumption that 
ordering of graphene planes within the electrode is pre- 
served over arbitrarily large distances. This assumption 
does not significantly affect any results, as we demon- 
strate in Sec. [Vj 

When the positive and negative electrodes are con- 
nected to opposite terminals of a voltage source, the ap- 
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FIG. 2. A schematic picture of a portion of the GS electrode, shown in a side view. The stack of graphene sheets (black lines), 
arranged with a separation c between layers, has some concentration of "spacer" molecules (tan-colored circles) in it. When a 
voltage is applied, ions (red circles with -|-'s) enter the electrode collectively in the form of disks whose size is determined by 
the average distance R between spacers in a given plane. 



plied voltage Vtot between them is related to the sys- 
tem's free energy Ftot according to Vtot — dFtot/dQ. In 
this paper, for simplicity, we focus only on the negative 
electrode, employing the definition of voltage given in Eq. 
([5]) . This voltage V is related to the total voltage applied 
between the electrodes according to Vtot = V + dF' /dQ, 
where F' is the free energy density of the configuration 
of Q/e anions and electron holes within the positive elec- 
trode. In this way the SC is equivalent to two capacitors 
connected in series, with the applied voltage Vtot being 
shared between the two halves of the system. One can 
also think that the results we derive below correspond to 
a device with a very large positive electrode, so that the 
smaller negative electrode's capacitance dominates the 
series sum and Vtot — V. 

When the voltage V is applied, cations are driven to in- 
tercalate into the negative electrode and they establish a 
uniform electrochemical potential throughout the device. 
In our treatment below we make the simplifying assump- 
tion that the central ionic liquid reservoir has a large 
negative chemical potential due to the strong. Coulomb- 
based correlation energy of ions within the ionic liquid 
reservoir, so that at V — both electrodes are empty of 
ions. In this case the capacitor acquires a finite charge 
only when V is larger than some positive "threshold volt- 
age" Vt- This behavior has been described previously for 
dual-graphite energy cells |29j . 

For practical SC devices, the effectiveness of the elec- 
trode is often limited by the strong binding of graphene 
sheets to each other through van der Waals (vdW) attrac- 
tion 131 [5] • The short-ranged vdW interaction between 
sheets produces a large contact energy 7 = (1.5 ± 0.4) 
eV/nm^ ^30j, and in the absence of any modification to 
the graphene this interaction results in dense, graphite- 
like agglomerates whose interior surface area is not read- 
ily accessible to ions. In this case only a large applied 
voltage can induce ions to enter the electrode by forcibly 
separating the space between adjacent graphene sheets. 
This process is often too slow to meet the fast charg- 
ing/discharging requirements of practical devices. Fur- 
ther, such agglomeration can require the device to oper- 
ate at large voltages that are outside the limited electro- 
chemical stability window of the electrolyte. 

One method of overcoming these difficulties is to fab- 



ricate electrodes from chemically modified graphene, in 
which graphene sheets are bonded to polymers or some 
other "spacer" molecule before being condensed into a 
stack similar to the one shown in Fig. [2j In this arrange- 
ment each spacer provides sufficient steric hindrance to 
separate adjacent graphene layers at a given point. Elec- 
trodes with such spacers have recently been constructed 
and shown to generate capacitance values in excess of 100 
F/cm'^ even at relatively low voltages 0. 

In the region surrounding each spacer, separated 
graphene sheets experience a strong elastic stress that re- 
sults from the competition between their attractive vdW 
interaction and the elastic energy associated with bend- 
ing them to bring them together. The introduction of 
intercalated ions into the GS requires additional separa- 
tion of graphene sheets, since the diameter a of the ions is 
also larger than the spacing c between adjacent graphene 
sheets (and is typically comparable to the size of the spac- 
ers). In equilibrium, this separation occurs in a way that 
best relieves the elastic stress surrounding the spacers — 
namely, through the formation of filled, 2D "pockets" of 
ions rather than through the introduction of scattered, 
individual ions. The size of these pockets is defined by 
the average distance R between spacers (see Fig. [2]). 

One can explain the formation of these pockets of ions 
using the language of elastic energy-mediated attraction. 
When a single ion enters the GS, it creates a region of 
intense elastic stress around itself. If two such ions are 
introduced to the GS, they can reduce the total elastic 
energy of the system by approaching each other closely, 
so that together they share the elastic energy associated 
with deforming the GS around them. Thus, there is an 
effective attraction between ions in the GS mediated by 
the elastic strain they induce. This attraction has been 
observed experimentally in intercalation experiments [191 
I3T] and has been described theoretically as an effective 
lateral attractive interaction [26l |27j . 

When there is some finite concentration of ions within 
the electrode, scattered ions between two adjacent sheets 
in the GS attract each other to form a large "disk" of 
ionic charge. The growth of the disk is truncated at 
the size i?, at which point ions fill the area between 
neighboring spacers. These spacers are generally mas- 
sive enough that they can be considered immobile on the 
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time scales of charging/discharging of the electrode. The 
2D charge density ctq of the disk is determined by the 
balance between the elastic energy-mediated attraction 
and the Coulomb repulsion between ions. 

The result is that, as a voltage is applied, the elec- 
trode is charged by the incremental addition of disks of 
ions (which may have a somewhat irregular shape) that 
enter the electrode collectively, as shown in Fig. [2j In our 
derivation of the capacitance in Sec. |V] we assume that 
the disk size R ^ a > c, so that every disk contains a 
large number of ions. 

One may well notice that in the limit where there are 
no spacers at all in the GS, the electrode is simply tur- 
bostratic graphite. In this case the disk size R oo, 
which suggests that ions enter the GS as infinite, uni- 
form planes. In fact, this phenomenon is well-known: it 
is referred to as "staging" of graphite intercalation com- 
pounds and has been studied for more than eighty years 
[m US]. In such compounds, guest ions such as Li+, 
K"*", or Ca^^ are forcibly introduced into pure graphite 
by an applied voltage. In most cases, the elastic energy- 
mediated attraction causes the intercalated ions to ar- 
range themselves in an ordered sequence of filled and 
empty interlayer galleries, with each filled ion layer hav- 
ing some fixed 2D concentration of ions and the distance 
between filled layers decreasing with the overall concen- 
tration of intercalated ions. The number of graphene 
sheets between two subsequent filled layers is called the 
"stage" of the compound: low densities of intercalated 
ions corresponds to large stage, while the maximum fill- 
ing of ions is called "stage 1." The concentration cro/e of 
ions within a filled plane depends in general on the ion 
size: small ions tend to have stoichiometry XCg in their 
stage 1 form, while larger ions form less dense arrange- 
ments such as XCg or XC12 [SI]. In each case, (To/(e/c^) 
is of order unity. The physics of staging and its implica- 



tions for capacitance are discussed more fully in Sec. IV 



is the quantum kinetic energy per unit area of electrons 
within the graphene and is the electrostatic energy. 

The quantum kinetic energy 'Wq can be obtained by 
integrating Eq. This produces (at T = 0) 



where 



(7) 



(8) 



is the chemical potential of graphene with the 2D elec- 
tron concentration rig = a /e. In terms of the capacitor's 
charge per unit area cr, becomes 



where 



(9) 



(10) 



is the dimensionless interaction constant of graphene. For 
free-standing graphene a « 2.2, so that for = 3.0 we 
have a ~ 0.7. 

Using the thermodynamic equations relating the en- 
ergy to the voltage V and capacitance per unit area 
produces an expression for the quantum capacitance per 
unit area "^g. In 2D the relations analogous to Eqs. ([s]) 
and ^ aieV = df^/da and = (d'^'W /da^)-\ so that 
the quantum capacitance is given by 



327ra2(£o£)V/e. 



(11) 



This quantum capacitance is added in series with the 
geometric capacitance 



<^,(d) 



(12) 



III. MEAN-FIELD THEORY OF CAPACITANCE 
IN A GRAPHENE STACK 

As a first attempt at describing the capacitance of a 
GS, one can try to use mean-field theory. This approach 
ignores any correlations among ions, such as the cluster- 
ing described in Sec.|Tlj In the simplest mean-field model, 
discrete ions are replaced by a uniform charge density 
eN+ that fills the electrode volume. This primtive pic- 
ture allows one to derive a value for the mean-field quan- 
tum capacitance Cq, which in the mean- field approach 
imposes an upper limit on the observable C{V). 

Before calculating the volumetric capacitance in the 
mean-field approach, one can consider first the problem 
of a single graphene sheet gated by a parallel metal elec- 
trode that is separated by a distance d from the sheet (as 
in the experiments of Refs. [El [16]). Such a system can 
be described as a parallel-plate capacitor, and its total 
energy per unit area is equal to '^q + ^ei, where '^q 



that results from the capacitor's electrostatic energy, 



a^d 



(13) 



Eqs. ([9| and (13) imply that when the charge of the 
capacitor is small enough that a e/9Tra^d^ we have 
^ "^eh and therefore the quantum capacitance dom- 
inates the series sum = {'rf~^ + ~ 'rfq. Con- 
versely, at cr e/dna'^d'^ we have 'rf ~'rfg. 

For a stack of many graphene sheets with a uniform 
compensating ion charge, the mean-field picture suggests 
that the electrode consists of many such single graphene 
capacitors connected in parallel. If one assumes that the 
distance between graphene sheets is a constant c, then the 
volumetric capacitance satisfies C{V) = '^(F)/c. Substi- 
tuting Eq. (Ill gives 



Cq{V) ~ 8a' 



V -Vt £o£r 
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(14) 
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In principle, this quantum capacitance should be added 
in series with a constant geometric term Cg that results 
from the system's electrostatic energy. However, the 
quantum capacitance Cq rises to the level of the geometri- 
cal value only at very large ion densities -/V+ ~ (a^c^)^^, 
which are not physically realistic. One can therefore 
say that the mean-field approach predicts a capacitance 
C{V) ~ Cq{V) over the entire relevant range of voltage. 
Cq iy ) is plotted as the thin hne in Fig. Is] 

We note here that Eq. ( 29 1 and all suEsequent formu- 



las for the volumetric capacitance are normalized to the 
volume of the empty GS, where graphene sheets are sep- 
arated by a distance c. 
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in Appendix \K\ 

In order to determine the energy of the staged arrange- 
ment of ions, one can approximate the ion layers as a 
sequence of uniform planes with surface charge density 
CTo positioned at 



z = (to + l/2)/i. 



0,±1,±2, 



(15) 



where h is related to the average ion concentration 
by 



h = cro/(eiV+) , 



(16) 



In this description the discreteness of graphene sheets is 
also ignored, so that h can take on values that are not 
integer multiples of c; this assumption is discussed at the 
end of this section. 

Within the TF approximation, the electrostatic poten- 
tial 4){z) and the 3D electron density Nf.{z) = nf,{z)/c are 
related by e4>{z) = ^e{ne{z)) + const. With the help of 
Eq. (|8| and a suitable choice of the additive constant, the 
Poisson equation for (t>{z) can be written as 



r{z) 



^ (^) 



(17) 



By the Gauss law, (j)'{z) exhibits a discontinuity of 
cro/(2eoer) across each ionic plane. 

For the case of large /i, the solution for 0(z) can be 
approximated by [13 [55] 



FIG. 3. The capacitance as a function of voltage for the 
hypothetical case where a = 1 and the disk charge density 
(To = e/c?, plotted for different values of the disk size R. 
At R/a ~ 1, the elastic energy-mediated attraction between 
ions is eliminated and the capacitance follows the mean-field 
quantum capacitance given by Eq. (14 1 (thin, black line). At 
i? — >■ oo, ions enter the GS in stages and the capacitance fol- 
lows the result of Eq. ( |24[| (dashed curve). The behavior for 
finite R, derived in Sec. |V[ is shown schematically by the red 
and blue thick lines for two different values of R. Proper nu- 
merical coefficients are shown on the axes for the thin black 
and dashed curves, while for the thick curves the location 
and height of the peak capacitance is indicated only by an 
approximate scaling relation. 



IV. CAPACITANCE OF STAGED GRAPHITE 

The mean-field picture presented in the previous sec- 
tion is apparently incompatible with the phenomenon of 
staging observed in graphite intercalation compounds, as 
described in Sec. |llj In staged graphite intercalated ions 
do not fill the electrode volume uniformly, but instead 
form a sequence of filled and empty interlayer galleries, 
periodic in the stacking direction z. To compute the ca- 
pacitance of this structure one can use the Thomas- Fermi 
(TF) approximation pSl . The key findings of this ap- 
proach are summarized below, with further details given 



(^) 



(18) 



2^2 47reoer (Az + zqY ' 
where Az is the distance from the nearest ionic plane and 
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(19) 



is the characteristic thickness of the plane's screening at- 
mosphere. The corresponding 3D electron density Ne{z) 
near ionic planes is given by 



^e(^) 



9c 



Ana^iAz + zoY 



(20) 



For large ionic plane separation h ^ zq one can talk 
about an effective repulsion of the planes due to 
the overlap of their screening atmospheres. The elec- 
tron density perturbation caused by the overlap is the 
most significant near the midplanes, e.g. at z = 0, where 
the density of the electrons and their screening ability 
are strongly diminished. Therefore, the interaction en- 
ergy u{h) per plane per unit area can be estimated as 
u{h) ~ e(j){0)Ne{0)h oc This estimate is verified by 

a detailed calculation (cf. Appendix [A]) , which provides 
the numerical coefficient: 



u{h) 



Cl 



Cl 



1.16953. 



(21) 
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The total energy per unit volume, 



(22) 



and one needs a more careful method to calculate the 
capacitance, as outlined in Appendix [A} T he results of 
these calculations are discussed in Sec. IVII 



includes contributions from both the interplane repulsion 
and the self-energy of each plane, which is 
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per ion (cf. Appendix [A]) and which enters into the 
threshold voltage Vf. Using these results, we can derive 
the voltage V and the capacitance C of the staged com- 
pound near the intercalation threshold as a function of 
interplane distance h. For large h this is done by substi- 
tuting the formula iV+ — ao/eh for the ion concentration 
and Eq. ( 22 1 for the energy density into the thermody- 
namic relations ([s]) and ([6]). Next, eliminating h, we ob- 
tain 



C{V) 



We should note that in the derivation of Eq. (24 1 it 
is assumed that the graphite can be represented as a 
continuous, electron-filled medium. In a more realistic 
treatment, one could treat h not as a continuous vari- 
able, but as a distance whose value is restricted to be 
an integer multiple of c, known as the stage number s in 
graphite [THl [H] • In this case the ideal periodic structure 
discussed above would exist only for a set of discrete ion 
concentrations = (Tq /{ecs). At such concentrations 
V would change discontinuously, giving a sharp peak in 
C{V). In between the peaks, the system would be a 
mixture of two stages with neighboring stage numbers. 
Equation (24) would then represent C{V) averaged over 



an interval of voltage containing several peaks. 

One could also notice that we have assumed ionic 
planes to be immersed in the continuum electron back- 
ground without any gap between the two. If nonzero 
gaps of width do ~ a/2 + c exist at both sides of every 
plane, then the energy density acquires an extra term 
eaoN-f-do/{eoSr)- Since (Tq is a constant, this extra term 
can be absorbed into Vfi 



Vt^Vt 



(25) 



Notably, the "geometric capacitance" [Eq. ( 13 )] only 
changes the threshold voltage but does not reduce C, 
unlike in the mean-field theory of Eq. ^ and Sec. 
The behavior of C{V) predicted by Eq. ([24| 
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markedly different from the mean-field theory Eq, 
as shown by the dashed line in Fig. [3j Most notably, the 
capacitance of staged graphite diverges at the threshold 
voltage rather than collapsing to zero, and it remains 
parametrically larger than the result given by Eq. ( 14 ) 
up until V — Vf ^ e/(47r£o£rc). At those large voltages, 
the approximate formula of Eq. ( 24 1 becomes inaccurate 



V. CAPACITANCE OF A GRAPHENE STACK 

For the main problem of this paper, as defined in 
Sec. [n) the different results of Eq. ^ and Eq. ([24]) 
present something of a puzzle. One the one hand, when 
the spacers in the GS electrode are plentiful, the graphene 
sheets are fully separated and ions are free to enter the 
GS individually and fill its volume uniformly. In this 
case, the capacitance should be similar to the mean-field 
result of Eq. (14). On the other hand, when the spac- 



ers in the GS are sparse, the intercalating ionic charges 
arrange themselves in the staged structure discussed in 
Sec. |IV| an d the capacitance should follow the prediction 
of Eq. (|24[). In this section we show that both relations 
C oc (V - Vt) [as in Eq. (Ill)] and C oc {V - 14)^*^^ N 



in Eq. (24)] are realized in suitable ranges of voltage and 



we explain how the crossover between them occurs. 

For the purpose of our conceptual discussion in this 
section we drop all numerical coefficients and focus in- 
stead on the scaling behavior of the capacitance as a func- 
tion of voltage. For simplicity, we consider the case where 
a ~ 1 and a/c is also of order unity. This second assump- 
tion amounts to assuming that cto ^ (^/a^) ~ (c/c^). 
The effect of allowing a or aQ/{e/c'^) to be small is dis- 
cussed briefly at the end of the section, but we note here 
that as long as a/c <C (R/c)^^'^ the scaling relations we 
derive below are unaffected. A more thorough derivation 
of the capacitance is given in Ref. [53t which includes an 
analysis of additional regimes that appear if a ^ 1 or 
Cto < e/c^. 

Our essential observation is that when disks of charge 
enter the GS, each disk draws around itself a strongly- 
bound, nonlinear screening atmosphere of electrons that 
renormalizes its charge. To see how this happens, con- 
sider first the case where a single disk is introduced into 
the GS. At very small distances z from the surface of the 
disk, the disk behaves like an infinite plane and the elec- 
tric potential can be described by Eq. (18). This formula 



for the potential is derived under the TF approximation, 
which neglects any quantum effects associated with fi- 
nite electron wavelength, and therefore its applicability 
is limited to regions where the Fermi wavelength Ai? is 
much smaller than the length scale over which the po- 
tential is varying in the transverse direction. Thus, to 
find the distance ztf at which TF screening ends, one 
can equate Af ~ 1/%/^ ^ e/(eo£r0) to the disk size R, 



which by Eq. ( 18 ) gives 



ZTF 



/Rc. 



(26) 



Since ztf ^ Ri the region over which the TF approx- 
imation is applicable can be described as a disk-shaped 
volume with thickness ztf ^ R and diameter R. 
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The net effect of the strong, TF screening at \z\ < ztf 
is to renormahze the charge of the disk. To see the ex- 
tent of this renormahzation, one can calculate the charge 
enclosed within the region \z\ < ztf- In this case the 
Gauss law can be written as q ^ —SQSrR^ {dcj) / dz)\z=z'j:-F , 
which gives q ^ e^jR/c. This charge q enclosed within 
the TF region is much smaller than the bare disk charge 
aoR^ ~ e{R/cf. 

Outside of the TF region, the potential created by the 
reduced charge q is not strongly screened. Indeed, at 
these large distances the potential is affected only by the 
relatively weak dielectric response of graphene layers with 
a small TDOS. Such linear screening is described in detail 
in Appendix |B] For the purpose of the present calcula- 
tion of the capacitance, however, it is sufficient to think 
that outside of the TF region there is no screening of the 
disk charge. Instead, the remaining negative charge —q 
required to make the system electroneutral is distributed 
uniformly throughout the electrode. 



linear 



-TF 



TF 



\ ^0 ^ 



R 



charge inside: q ~ e (R/cY'^ 



FIG. 4. A schematic depiction of the screening atmosphere 
surrounding a disk of ionic charge (rectangle) with diame- 
ter R inside the graphene stack. The limit of TF screening, 
which corresponds to a distance ztf from the disk, is shown 
by the dashed hne. The region outside the TF region corre- 
sponds to linear screening. The characteristic decay length 
zo of the nonlinear potential is indicated by the dotted line. 
The total amount of charge inside the TF screening atmo- 
sphere is denoted q and is much smaller than the bare charge 
ao-R^ ~ e{R/cf' of the disk of ions. 

We can now consider what happens when there is some 
finite concentration N of these charge-renormalized disks 
within the volume of the electrode. When the concentra- 
tion of the disks is low enough that N <^ 1/{R-^ztf), the 
disks are well-separated from each other and their non- 
linear screening atmospheres do not overlap. One can 
think, then, that in this low density limit the capacitor 
charge consists of a dilute gas of disks, each with effec- 
tive charge q, neutralized by a uniform background of 
electronic charge with charge density —qN. 

In this configuration, the total capacitor energy U has 
three components: the self-energy of the disks with ef- 
fective charge q, the electrostatic energy Uei associated 
with the Coulomb interaction between the disks and the 
uniform background, and the quantum kinetic energy Uq 



of the uniform, neutralizing electronic charge. The self- 
energy component affects only the threshold voltage Vt 
and does not enter into the capacitance. The Coulomb 
energy Uei is much smaller in magnitude than Uq in the 
limit N ^ l/(i?^ZTF); this is shown explicitly in Ap- 
pendix |B] Thus, in the regime where the disk screen- 
ing atmospheres do not overlap, the capacitance is domi- 
nated by the quantum kinetic energy of the uniform elec- 
tronic background. 

The quantum kinetic energy Uq can be derived in a 
way that is similar to the derivation of Eq. ^ . Here, the 
3D electron density (in regions outside the TF screening 
atmospheres of the disks) is given by = qN/e. Us- 
ing q ^ e^/R/c and Eq. ([7| for the energy of graphene 
per unit area gives for the energy per unit volume Uq ~ 
(7Vc3)3/2(i?/c)3/4(e2/£o£^c4). Since is related to the 
total capacitor charge Q per unit volume hy Q = aoR^N, 
one can use the thermodynamic relations for capacitance 
and voltage in Eqs. ([s]) and (|6| to get 



C{V) 



9/2 



V -Vt eoEr 



e/ilTSQSrC 



(27) 



if -Ta — ^ 



Corrections to this result associated with the Coulomb 
interaction between disks are discussed in Appendix |B] 
[Eq. ~ 



In Eq. (27) one can see the strong effect of the renor- 



mahzation of the disk charge. This equation is similar 
the mean-field relation of Eq. ( 14 1 in the sense that both 
give C oc {V — Vt), but Eq. (27) has a significantly larger 
coefficient. Indeed, the slope of the C{V) relation in 
Eq. ( p7| is larger than that of the mean-field result by 
a factor {R/c)^/"^ ^ 1. This large enhancement of the 
capacitance can be viewed as a direct result of the non- 
linear screening of each disk. The strongly-bound TF 
screening atmosphere surrounding each disk greatly re- 
duces the uniform electron concentration in the electrode 
for a given capacitor charge Q and therefore leads to a 
smaller quantum kinetic energy cost for capacitor charg- 
ing. 

When the capacitor charge is made larger, such that 
the concentration of disks A^ ^ l/(i?^ZTF): the nonlin- 
ear screening atmospheres of adjacent disks overlap and 
one can no longer talk about a uniform electron charge 
—qN filling the space between disks. Instead, the capac- 
itance is dominated by the repulsive interaction between 
neighboring disks, and the capacitance is described by 
the staging theory of Eq. ( 24 1 : 



C{V) 



e/47reo£rc\ "^^^ eo£r 
V-Vt 



if yyCL-^'^ 



(28) 



The resulting CiV) dependence, which combines 
Eqs. (27) and (28), is shown schematically in Fig. [sj No- 
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tice that in the derivation of these results it was not nec- 
essary to assume ordering of the graphene sheets over 
distances larger than R. Indeed, the low-density re- 
sult of Eqs. (27 1 assumes only that the concentration of 



graphene sheets is roughly uniform throughout the GS, 
while Eq. ( 28 1 is based on disks interacting over a dis- 
tance z <^ ztf <C R. Therefore, our earlier assumption 
of long-range ordering of the GS does not significantly 
alter any results. 



In this way the crossover between Eqs. (14 1 and (24) 



can be understood as follows. At low voltages, one can 
still think that the result C cx (V—Vt) is a consequence of 
uniformly raising the electron Fermi level throughout the 
GS in order to provide a neutralizing electron concentra- 
tion. However, this neutralizing concentration should be 
thought of as a compensation not to the total ionic charge 
but to the much smaller renormalized ionic charge. It is 
this renormalization that allows the capacitance to be 
large and produces a smooth crossover to the behavior 
C oc — Vt)~^^^ associated with staging. 

The most dramatic consequence of this renormaliza- 
tion is the large peak of the capacitance at low voltage 
{V -Vt)/{e/'iiTeoerc) ~ {c/Rf/^, as shown in Fig.js] In- 
deed at this point the capacitance attains a value which 
is larger than the maximum mean-field capacitance by 
the parametric factor (i?/c)^. This factor can explain 
the large difference between the mean-field estimate and 
the large observed values of volumetric capacitance. 

In deriving the capacitance maximum we assumed that 
a and ao/{e/c^) ^ {c/aY were both of order unity. If a 
or CTo are reduced, as should be expected when a > c, 
then naturally this peak in capacitance declines, since 
reducing either a or ctq implies a weaker role of Coulomb 
interactions relative to the quantum kinetic energy and 
therefore a weaker renormalization of the ion charge. 
In fact, one can show |33] that as long as a {c/R) 
and (To/(e/c^) ^ (c/R)^^^, the peak in capacitance of 
Fig. ^should be replaced by the somewhat smaller value 
~ a'^agc^ /e)'^{R/c)'^{eoer/c'^), while the overall qualita- 
tive picture of C{V) is not affected. 

Finally, one can notice from Fig. [3] that a sparser ar- 
rangement of spacers within the GS, which corresponds 
to larger R, results in larger capacitance and therefore 
in greater energy storage for a given voltage. However, 
such increased capacitance comes at the cost of slower 
capacitor charging. That is, when spacers within the GS 
are sparse the process of ion intercalation into the GS 
is slow kinetically, as discussed in Sec. |TT] and therefore 
the power of the device is reduced. One can thus say 
that there is a fundamental tradeoff between high energy 
density and high power in graphene supercapacitors that 
can be adjusted by altering the density of spacers within 
the GS. 



VI. DISCUSSION 

In the previous section we showed how the capacitance 
can be much larger than the mean-field quantum capaci- 
tance Cq as a result of the elastic energy-mediated attrac- 
tion between ions and the nonlinear screening of disk-like 
ion bunches by the surrounding graphene layers. These 
effects together produce a large, i?-dependent peak in the 
capacitance at small voltages, as shown schematically in 
Fig. HI 

At larger voltages the capacitance is determined by 
the physics of staged graphite, as described in Sec. |IV| 
In the present section we discuss in greater detail the ca- 
pacitance at these large voltages, and we show that near 
the steric limit of capacitor charging, where neighbor- 
ing disks are separated by only a few graphene sheets, 
there is a noteworthy deviation of C{V) from the result 
of Eq. ( [24| ) toward larger capacitance. This deviation 
brings the schematic picture of Fig. [3] closer in line with 
the large, mostly-constant C{V) curve observed experi- 
mentally. 

In our derivation of the approximate C{V) relation 
of Eq. (24), it was assumed that disks of ionic charge 



are separated by a distance h ^ zq, so that their TF 
screening atmospheres overlap only weakly. This as- 
sumption becomes invalid when the capacitor charge is 
large enough that neighboring disks are separated by only 
a few graphene layers. Instead, to calculate the capaci- 
tance at such large voltages one should employ the full 
solution of the TF equation [Eq. (ITt])] rather than the 



asymptotic interaction law given in Eq. (21 1. Such 
a calculation is presented in Appendix |A| and the re- 
sults for capacitance and the stage number s are shown 
in Fig. [5] for = 3 and cto — e/(lnm)^. In this plot 
the capacitance is normalized to the volume of the filled 
electrode at stage s = 1 (which includes the volume of 
the ions as well as the GS itself), and assumes an ion 
size a — 1 nm « 3c. The right-most point of Fig. [5] cor- 
responds to stage 1, after which presumably no further 
capacitor charging is possible. 

One can notice from Fig. [5] that at small V ~ Vt the 
capacitance declines with voltage accordin g to the depe- 
dence C oc {V - Vt)~'^^^ derived in Sec. |iv) At larger 
voltages, however, the capacitance attains a minimum 
and then increases weakly with voltage. This deviation 
from Eq. ( 24 ) allows the capacitance to be larger and rela- 



tively constant over the majority of the operating voltage 
range. 

The origin of the upward deviation from Eq. (24) 
can be understood as follows. At large enough volt- 
ages that /i < zq, the disks' screening atmospheres over- 
lap strongly and the electron density becomes roughly 
uniform throughout the GS. As a result, the energy at 
these large voltages is dominated by the quantum kinetic 
energy associated with uniformly raising the graphene 
Fermi level. The capacitance therefore behaves simi- 
larly to the mean-field quantum capacitance described 



by Eq. (14), rising linearly with voltage. The linear rela- 
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FIG. 5. (Color online) The volumetric capacitance C (thick 
black curve, left axis) and the stage number s (thin blue 
curve, right axis) as a function of voltage in the Thomas- 
Fermi approximation for turbostratic graphite {R = oo), with 
(70 = e/(lnm)^, and e,, — 3. The capacitance in this plot is 
normalized to the volume of stage- 1 graphite with interca- 
lated ions of diameter a — 1 nm. The dashed black line is an 
approximation for C{V) at large voltage, given by Eq. ||29[). 



tion describing this portion of the C{V) curve is shown 
as the dashed line in Fig. [5] and is derived below. Its 
relatively small slope is identical to that of the mean- 
field quantum capacitance equation in Eq. ( |14[ ), since at 
these large ion densities there is no renormalization of 
the ion charge. Crucially, however, there is effectively a 
positive vertical offset of the quantum capacitance curve 
that allows the capacitance to be large. 

The value of this vertical offset can be calculated by 



first noting that the addition of a constant to Eq. ( 14 ) is 



equivalent to substituting a shifted value of the threshold 
voltage Vt ■ This shift in Vt can be understood intuitively 
as the effective elimination of the self-energy of each disk. 
That is, at /i < zq each disk no longer has a well-defined 
screening atmosphere, so that there is no concept of a 
constant self-energy contributing to the threshold volt- 
age. Thus, an equation for the dashed line of Fig. [5] can 
be obtained by replacing Vt with Vt — uq/b in Eg. (|14| , 
where uq is the self-energy per ion of an isolated plane 
in the GS [Eq. (23)]. Making this substitution gives the 
following expression: 



C{V) 



/207367r2 



1/3 



8a' 



125 
V-Vt 



£o£r 



2/3 



(29) 



For situations where zq is as large as a few times c, such 
as may result when ions are large enough that ctq is rel- 
atively small, Eq. (29) can occupy the majority of the 
SC's operating range of voltage. This may help to ex- 
plain why experiments on graphene SCs, where z^jc is 
probably between 1 and 3, generally report a capacitance 
that varies only slightly with voltage and retains a value 
of ~ 100 Y jcvc? over the entire operating range of voltage 
[3Hi[28j. 



At very small voltages, on the other hand, one should 
still expect a strongly-increasing, linear C(T^) relation as- 
sociated with the quantum capacitance of the renormal- 
ized disk charge, as discussed in Sec. |V] This behavior 
leads to a large peak in the capacitance at a particular 
small voltage. The apparent lack of such a sharp peak in 
experimental data is likely the result of finite tempera- 
ture or disorder in the graphene stack. Both disorder and 
finite thermal energy of ions work to diminish the posi- 
tional correlations among disks of ionic charge when these 
disks are distant from each other, leading to a larger av- 
erage interaction energy between disks and therefore to a 
smaller, somewhat smeared capacitance peak. Nonethe- 
less, even if positional correlations between disks are lost 
at stage s > 4 or so, as is commonly reported for tradi- 
tional graphite intercalation compounds [31] , Fig. [5] sug- 
gests that a low- voltage capacitance peak of several hun- 
dred F/cm^ may still be observable. 

Long-range elastic interactions between ion bunches 
may also play an important role for the capacitance, and 
these will be explored in a later publication. 



ACKNOWLEDGMENTS 

This work was supported primarily by the MR- 
SEC Program of the National Science Foundation un- 
der Award Number DMR-0819885. Additional support 
comes from NSF Grant No. PHYG5-51164 and from 
UCOP. 



Appendix A: Thomas-Fermi theory of the staging 
regime 

In this Appendix we provide the details of the TF the- 
ory of intercalated graphite discussed in Sec. |IV| It is 
convenient for our treatment here to introduce dimen- 
sionless variables (/>(z) = e0(z)/i?o, z — zj zq = z^jt^ 
and h = h/ where 



(Al) 



are the energy and length units, and 

1/6 



/ 



97r (To 2 2 

a c 

4 e 



(A2) 



is a dimcnsionless parameter. In further calculations, it is 
sufficient to consider a single period —h/2 < z < h/2, so 
that Eq. (17) yields the following boundary- value prob- 
lem: 



(A3) 
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Its solution can we represented by the inverse function 



which is equivalent to Eq. (23). The exact numerical 



I{x) 



du 



I) 



(A4) 

where (j)o stands for 0(0). Note that I{x) can be ex- 
pressed in terms of the Gauss hypergeometric function 
2Fi{a,b; c; x) and the Euler gamma function r(a::): 



coefficient in the interplane interaction energy u{h) per 
unit area [Eq. (21)] is Ci ~ I^/{60tt^), which is about 
10% smaller than what was obtained in Ref. 

In order to compute the voltage V we take the deriva- 
tive of F with respect to the ion concentration: 



eV 



dF _ dF/d(j)o 
dN+ ~ dN+/d^o 



(A13) 



' 6' 2' 6' a;3 



'6n 



r(7/6) 
r(2/3) 



2.97448 . 



(A5) 
(A6) 



This concentration can be expressed in terms of our vari- 



ables using Eq. (A4). The result is 



For distant planes, /i 3> 1, we have (/)/0o 3> 1. Expanding 
I{(f>/<f>o) in Eq. (A4) to the leading order in this ratio and 



^ eh' c ''y2Pcbo 



(AM) 



going through a simple algebra, we can show that near 
the ionic planes the potential behaves as 



where s is the dimensionless "stage number." Equa- 
tions (|A11|), (|A13|), and iAUl imply 
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(A7) 



V = 




(A15) 



which is equivalent to Eq. ( |A7[ ). 
middle point z = is given by 

00 ^ ^llolh" 



The potential at the Taking another derivative, we get the capacitance: 



(A8) 



which implies that at this point Eq. ( A7 ) errs by the fac- 
tor /^ / 6 « 1.47. This deviation is caused by the overlap 
of the screening atmospheres of the planes mentioned in 
Sec. |IV| The 3D electron density in physical units can be 
computed from 



iVe(z) 



02 (^) 



(A9) 



which leads to Eq. ( 20 ) 



Let us now compute the capacitance. First, we need to 
calculate the free energy density F of the system. Keep- 
ing the electrostatic energy and electron kinetic energy 
but neglecting the entropy, we have 



dN+ ,aoIhcbl + 

Cq = e —— = 27reoerQ! — ■ 
dV e 



(A16) 



It is easy to check that in the limiting cases of high and 
low ion concentrations we recover Eqs. (11) and (24 1, 
respectively. 



By virtue of Eqs. (|A5j), (|A6|), {Ku\ , and ( |AT4| -( |A16| , 
all quantities of interest are functions of 0o- It is then 
possible to graph the dependencies of the capacitance Cq 
and the stage number s on voltage as parametric plots. 
An example is shown in Fig. [5] and discussed in Sec. |VI[ 



Appendix B: Linear screening and renormalized 
mean-field theory in a graphene staclc 
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Evaluating the integral following Ref. | 



F 



0/1 



we get 



bee \h 3 

(All) 

where 0/i denotes the potential at the ionic plane: 0/i = 
(j){h/2). Below we also use the shorthand notation = 
I {iph / 4'o) ■ Considering again the limit of large interplane 
distance, we can write the result of Eq. ( All ) as Eq. ( 22 ) 
with 



In this Appendix we compute the screened potential of 
a Coulomb charge in the undoped GS and use it to calcu- 
late the correction to the renormalized mean-field theory 
expression for the volumetric capacitance [Eq. (27)]. We 



assume a relatively sparse filling of the GS by ions, so 
that the volume of the GS is not significantly expanded 
and we can still think of the GS as a stack of graphene 
sheets with separation c. 

Due to the anisotropy of the system the screened elec- 
tric potential 0s(p, z) surrounding a point charge q is 
anisotropic as well. Here p is the radial coordinate in 
the x-y plane. The starting point of the calculation is 
the electron dielectric function of the GS, which can be 
written as 



eVt 



Eof, 



(A12) 



,ik,k,) = l-e(f>c{k,h)P{k,k,), 



(Bl) 
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where k ~ ^/fc^ + is the in-plane momentum, 



0c (fc, kz) 



SoSr 



1,2 



(B2) 



is the Coulomb potential in the Fourier space, P{k,kz) 
is the polarization function, and henceforth the tilde 
marks the Fourier transform of the corresponding quani- 
tity without the tilde. 

If adjacent graphene layers are assumed to have no 
electronic interlayer coupling, as in turbostratic graphite, 
then P{k, k^) is related to the polarization function of 2D 
monolayer graphene, P2, simply by P{k,kz) = P2{k)/c. 
The polarization function P2{k) has been previously cal- 
culated BSl ESI to be 



4 hv 



(B3) 



With this result one can define the potential (j)s{k,kz) in 
Fourier space: 



0s (fc, k^ 



(t)cik,k^) _ q 

ee{k,kz) eoErlk'^ + k'^ + na\k\/c] 

q 

£o£r[fcz + 7ra|fc|/c] 



(B4) 



The last, approximate relation in Eq. (B4) is valid for 
distances p c, so that k^ <^ \k\/c. Taking the Fourier 
transform of this equation gives the following result for 
the potential: 



(f>s{p, z) 



47reo£rC 



A-2/3P I 4-2/3: 
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(B5) 



where Ap and A 



given by Ap 



are numeric coefhcients of order unity, 
r(3/4)V%/^ and A, = A/{TT^a^). 



For the case of a disk of charge within the GS, Eq. (B5 ) 



can be used to describe the potential outside the TF 
screening atmosphere by substituting for q the renormal- 
ized disk charge q ^ e^J R/ac. It is worth noting that, to 
within numerical coefficients, Eq. (B5) smoothly matches 



gion 
p = 



the TF result of Eq.JlSf at the boundary of the TF re- 
That is, Eqs. (18 1 and (B5| are equal at the points 
z — zxF ^ \J Rcja and z = 0, p ^ i?, which lie 
on the same equipotential contour. Thus the potential 
surrounding a charged disk can be described by the TF 
result of Eq. (18 1 at \z\ < ztf, p < R and by the Hnear 



response result of Eq. (B5) otherwise, with no parametric 
intermediate regime. 



associated with the configuration of positively-charged 
disks residing on a negatively-charged background. This 
approach was justified because when the concentration of 
disks N is smaller than 1/(_R^ztf) the Coulomb interac- 
tion energy between disks is much smaller than the quan- 
tum kinetic energy associated with the uniform electron 
charge. In the remainder of this Appendix we explicitly 
calculate the Coulomb energy and prove this inequality. 
We also find the small correction to the capacitance as- 
sociated with the disks' Coulomb interaction; this is pre- 
sented in Eq. (|B8 I . 



As explained in Sec.|Vj when the concentration of disks 
is very small the capacitor charge consists of a sparse 
arrangement of charge-renormalized disks with charge q 
and concentration N surrounded by a uniform electron 
charge with density —qN. The quantum kinetic energy 
per unit volume associated with the uniform electron 
charge is roughly 



R 



3/4 



(B6) 



as presented in Sec.jVj In the remainder of this Appendix, 
as in Sec. [Vj we drop all numeric coefficients and focus 
instead on parametric dependencies. We also again as- 
sume that a ^ 1 and ctq ~ e/c^; the general case of small 
a and ctq is examined in Ref. 33; 

In order to estimate the magnitude of the Coulomb en- 
ergy, one can consider that in their lowest energy config- 
uration, the disks form a correlated arrangement on the 
uniform background such that the disks minimize their 
repulsive energy while maintaining the fixed concentra- 
tion N . This arrangement is characterized by the average 
spacing between disks in the p and z directions, which we 
denote dp and dz, respectively. The repulsive interaction 
between disks is dictated by the linear potential given in 
Eq. (B5|. Since this potential is anisotropic, one can ex- 
pect that the spacing between disks is also anisotropic. In 
other words, the minimum energy arrangement of disks 
is that of an anisotropic Wigner crystal. (This situation 
is similar to the better-studied system of colloidal parti- 
cles that form a charge-renormalized 3D Wigner crystal 

The distances dp and dz can be found by noting that, in 
their minimum energy configuration, disks are arranged 
within the GS so that all nearest-neighbor interaction 
energies are equal in magnitude. This implies that dp and 
dz are determined by the relation (f>s(dp,0) ^ (j)s{Q,dz). 
Since the concentration of disks N ^ (d'j^dz)^^, one can 
solve for dp and dz as a function of N. This process 
gives 



dp/c ~ (7Vc3)-2/5 and dz/c - (Nc^)-^/^, so that 



Using the linear potential of Eq. ( B5 ) , one can calcu- the typical nearest-neighbor interaction energy is u„ 



late the Coulomb energy associated with a finite concen- 
tration of disks that are sufficiently separated from each 
other that their TF screening atmospheres do not over- 
lap. This Coulomb energy was ignored in our calculations 
of the capacitance in Sec.|V]at small V — Vt- Indeed, the 
expression of Eq. (27) is based on the quantum kinetic 



energy of electrons and neglects the electrostatic energy 



#,(dp,0) = qcj)Mdz) ~ (A^c3)3/5(i?/c)(eV£o£c). 

From the nearest-neighbor interaction energy one 
can estimate the total Coulomb energy of the anisotropic 
Wigner crystal. This Coulomb energy is, in fact, nega- 
tive, as in the case of an ordinary isotropic Wigner crystal 
[55] . since the attraction of each disk to its Wigner-Seitz 
cell of negative background charge is stronger than the 
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repulsion between neighboring disks. The magnitude of 
the Coulomb energy per disk is of the order of the nearest- 
neighbor interaction energy Mhu, and therefore the total 
Coulomb energy per unit volume Uei ^ —Nu^^, so that 



MS/5 



(B7) 



Comparing Eqs. (B6) and (B7) suggests that \Uei\ ^ 
\Uq\ whenever N < yJl/cR^ ~ II{R^ztf)- Thus, our 
assumption that the Coulomb interaction is unimportant 
for the main term of the capacitance is justified. 



Since the Coulomb energy of Eq. ( B7 1 is parametrically 



smaller than the quantum kinetic energy over the rele- 
vant range of voltage, {V - Vt) / {e/ 'iTreoSrc) <C (c/R)^^'^, 
its affect is only to provide a small correction to the main 
term of the capacitance, Cq, which is given by Eq. (27). 



Specifically, since the Coulomb energy is negative and 
small, one can say that it produces a large negative ca- 
pacitance per unit volume that is added in series with the 
relatively smaller main term. This negative capacitance 
is an extension of the well-known negative compressibil- 
ity of a conventional Wigner crystal 05] , and it produces 



a small positive correction to Eq. (27 1. Taking the ap- 



propriate derivatives of the total energy Uq + Uei gives 



C^Ca 



5/2 



V~Vt 



1/5- 



(B8) 



where Cq is given by Eq. (|27|) . Notice that the correction 



term in Eq. ( B8 ) grows to order unity precisely at the 
crossover point (V — Vt)/(e/47reo£c) ~ (c/R)^^^, where 
the TF screening atmospheres of neighboring disks begin 
to overlap and the capacitance transitions to the "stag- 
ing" resuh C cx (y - Vt)-^/^ of Eq. 



It should be emphasized that the above derivation of 
Eq. (B8) is schematic, and misses any numeric coeffi- 



cients multiplying the Coulomb correction. Such numeric 
coefficients may increase the magnitude of the Coulomb 
correction at small voltage and are potentially quite im- 
portant. A more careful calculation is the subject of a 
later publication. Nonetheless, our main conclusion that 
C should collapse to zero at {V — Vt) ~ 0, as does the 
quantum capacitance Cq, remains valid. 
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